Self-trapping in an array of coupled ID Bose gases 
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We study the transverse expansion of arrays of ultracold 8 'Rb atoms weakly confined in tubes created by 
a 2D optical lattice, and observe that transverse expansion is delayed because of mutual atom interactions. 
A mean-field model of a coupled array shows that atoms become localized within a roughly square fort-like 
self-trapping barrier with time-evolving edges. But the observed dynamics is poorly described by the mean- 
field model. Theoretical introduction of random phase fluctuations among tubes improves the agreement with 
experiment, but does not correctly predict the density at which the atoms start to expand with larger lattice 
depths. Our results suggest a new type of self-trapping, where quantum correlations suppress tunneling even 
when there are no density gradients. 



Ultra-cold atomic gases trapped by light are well- 
characterized many-body quantum systems. Without the dis- 
order that is common in condensed matter systems, theoretical 
analyses of cold gases in equilibrium can be extremely accu- 
rate, although high powered numerical techniques must some- 
times be employed Q). Cold gas experiments are also well- 
suited to studying out-of-equilibrium dynamics, like the evo- 
lution of many-body correlations, since experimental time- 
scales are relatively slow yet faster than typical relaxation 
and decoherence rates 2^4]. Significant progress in the de- 
scription of the out-of-equilibrium dynamics of ID bosonic 
and fermionic systems has been achieved thanks to the exis- 
tence of exact solutions [5 1 and the development of numerical 
methods such as time dependent density-matrix renormaliza- 
tion group (t-DMRG) |6| and time-evolving block decimation 
(TEBD) [7 1 methods. Those approaches fail in 2D and higher 
dimensions, where most out-of-equilibrium calculations rely 
on mean field, phase space methods, or approximate analyti- 
cal techniques which are generally restricted to the weak in- 
teraction regime [8 |. The experiment-theory comparisons pre- 
sented in this paper illustrate the breakdown of these approxi- 
mate approaches, highlighting the need for more computation- 
ally tractable methods for dealing with intermediate coupling. 
The particular way the theory deviates from the experiment 
strongly suggests a qualitatively new non-equilibrium effect. 

Macroscopic self-trapping in quantum degenerate Bose 
gases, where mean-field energy gradients suppress tunnel- 
ing, presents an interesting set of non-equilibrium phenom- 
ena, with analogs in non-linear photon optics |9l [Toll and 
Josephson junction arrays HI IB . Self-trapping has been stud- 
ied theoretically in all dimensions [ 1 2—20 1 and experimentally 
in double-well systems ll2"Tl and in arrays of 2D pancake-like 
Bose-Einstein condensates (BECs) created by a deep ID lat- 
tice potential ll22l . These self-trapping experiments are in the 
weak coupling limit, where mean-field theory clearly applies. 
In this letter, we experimentally and theoretically investigate 



self-trapping behavior in an array of coupled quasi- ID tubes 
created by a 2D optical lattice. Atoms freely expand along the 
axis of the tubes, so that their density decreases with time, un- 
til they eventually become too dilute for self-trapping and ex- 
pand ballistically transverse to the tubes. In contrast to previ- 
ous self-trapping work, our quasi- ID gases are in the interme- 
diate coupling regime ll23ll . To provide a baseline theoretical 
description of our experiments, we build a mean-field model 
of expanding coupled ID gases. We then incorporate fluctu- 
ations into our mean-field treatment via an approximation to 
the so called Truncated Wigner Approximation (TWA) 1S1 I241 
by introducing random tube-to-tube phase fluctuations with a 
tunable magnitude at the initial time of the evolution. At low 
lattice depths, small phase fluctuations improve the agreement 
with the experiment, while the dynamics in deeper lattices 
are better-described when the phases are maximally random- 
ized. Still, the self-trapping seen in the experiment is more 
widespread across the array and persists to lower densities 
than mean-field based models predict. 




line of sight t (ms) 



FIG. 1: (Color online) (a) Experimental setup (not to scale): Atom 
clouds are confined in a crossed-dipole trap and partitioned by a blue- 
detuned optical lattice into a coupled array of vertical tubes, (b) 
Peak transversely integrated axial density as a function of time for 
Vo = 13Er. The linear density drops rapidly with time as the atoms 
expand axially, mostly driven initially by their mutual interaction en- 
ergy. 
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It is instructive to first describe the pure mean field theory 
evolution qualitatively. When the BEC is loaded into the 2D 
lattice, the atoms are distributed among tubes with an approx- 
imately Thomas-Fermi profile. The tube-to-tube mean-field 
energy difference grows with the transverse density gradient. 
For a given peak density and tunneling energy, there is a crit- 
ical distance along the lattice directions at which the mean- 
field imbalance between adjacent tubes is too large for tun- 
neling to conserve energy. Atoms at larger radii will not tun- 
nel radially outward, while atoms at smaller radii will start to 
expand, at some point reflecting from the self-trapped edge. 
A similar thing happens in ID self-trapping [22], but there 
are more features in the 2D generalization. Atoms at 45° 
to the lattice axes are only self-trapped at larger radii, since 
their density gradients are smaller in the lattice directions. 
The distributions therefore develop a density depression in 
the middle with a roughly square fort-like barrier around the 
edges. Since all atoms tunnel in some lattice direction, the 
self-trapped edges evolve with time. In the absence of axial 
trapping, densities, and hence tube-to-tube density gradients, 
drop. Self-trapped edges are lost and the cloud ultimately ex- 
pands ballistically in the (x, y) plane. 

Details of the experimental setup are given in 11251 . and the 
experimental geometry is depicted in Fig. A BEC with a 
barely detectable impurity fraction and N w 3.5 x 10 5 87 Rb 
atoms in the \F = 1, mp = 1) state is produced in a crossed 
optical dipole trap, created by the intersection of two 1.06 /im 
wavelength horizontal laser beams. The trapping frequencies 
experienced by the BEC are u) x>y = 2ir x 38 Hz in the hori- 
zontal plane, and ui z = 2ir x 94 Hz vertically. A vertical mag- 
netic field gradient of 30.5 G/cm levitates the atoms. A two- 
dimensional square optical lattice is created using two slightly 
different frequency pairs of linearly polarized retro-reflected 
beams with 1 /e 2 widths of 650 and 705 /im, blue-detuned by 
5.2 THz from the D2 transition in 87 Rb. The 2D lattice is 
ramped up in time according to s(t) = 1/(1 + at 2 ), where 
a = 4.15 ms is the time constant. Non-adiabaticity in the 
lattice turn-on, measured by turning off the lattice with the re- 
verse ramp, adds less than 4 nK-fcs of energy, which is small 
compared to the ID chemical potential (150-250 nK-ks)- 

After the lattice is turned on, the crossed dipole trap beams 
are suddenly turned off so that the atomic distribution evolves 
in the 2D lattice alone. After a time t, we image the in-situ 
spatial distribution of the cloud using high-intensity absorp- 
tion imaging l26l . This technique involves illuminating the 
atoms with a resonant probe beam with intensity / 7 sat , 
where 7 sat is the saturation intensity. After the beam passes 
through the atomic cloud, it is blocked by a 400 /im diam- 
eter dark spot in the Fourier plane of a one-to-one imaging 
system. Using Babinet's principle, one can show that the de- 
tected signal is proportional to the square of the atomic den- 
sity distribution integrated along the line of sight |26|. For 
the highest densities there is some lensing of the absorption 
beam, which we correct for by observing the axial expansion 
in a 35Er lattice, where there is no transverse expansion. At 
our highest density, this correction is 10% of the measured 



root-mean-square (RMS) width. 
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FIG. 2: (Color online) Transverse RMS widths of the density dis- 
tributions as a function of time for (a) Vo = 7.25 Er, (b) 9.25 Er, 
(c) 11 Er, (d) 13 Er. Each data point is the average of 10 measure- 
ments. The error bars are the quadrature sum of the statistical noise 
and the uncertainty from the density dependence of the absorption 
beam lensing. There is an overall, systematic uncertainty associated 
with the imaging resolution of 0.5 /im; comparable shifts of the ini- 
tial size simply shift the theoretical curves at early times (below 10 
ms), followed by a gradual divergence of the curves at longer times. 
The blue dashed line is the result of a mean-field calculation of the 
dynamics with no random phase between the tubes (77 = 0). The red 
dotted line is for 77 = 0.2, the green dash-dotted lines are for rj = 0.4 
(a) and 77 = 0.5 (b), and the orange solid line is for 77 = 1. 

We have taken images of the evolving atom distributions 
from along the z-axis, where the self-trapping depressions 
would be directly observed. We do not see them, although 
calculations show that our imaging system could resolve them 
if they were there. Instead, when viewed from along z or 
transversely 45° from the lattice axes, the integrated density- 
squared profiles are featureless and well-fitted by Gaussians. 
All data displayed below is taken from the side view. We in- 
tegrate the distributions over the central ±10% of the vertical 
Thomas Fermi radius, so that the atoms we consider within 
each tube have approximately the same density. We fit the 
resulting ID distribution to a Gaussian and extract a n 2, the 
transverse RMS width. In Fig. [TJ? we plot an example of how 
the transversely-integrated axial peak density varies with time, 
which we derive from measuring the axial expansion of atoms. 
In Fig. [2] we display a n 2 as a function of evolution time in the 
2D lattice for lattice depths of V = (7.25, 9.25, 11, 13) E R 
(uncertainty ±2.5%), where Er — h 2 k 2 /2M is the recoil en- 
ergy and M is the atomic mass. 

The behavior of a 2 n is qualitatively different for Vq = 
7.25 Er compared to larger lattice depths. For Vb = 7.25 Er, 
the RMS width of the cloud increases quickly from t = and 
slows near t = 5 ms, where the curve changes its concav- 
ity. For V a = 9.25, 11, and 13 Er, the RMS width remains 
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constant to within experimental uncertainty for t < t c , after 
which it increases linearly in time. The time t c at which the 
cloud begins to expand increases as the lattice depth increases. 
As noted earlier, the density decreases with time as the atoms 
expand axially in the tubes. We can thus associate with the 
time t c a 3D density at the center of the cloud, p(t c ). We 
find that the ratio p{t c ) /J — {1190 ± 310, 860 ± 190, 910 ± 
210} pm^Ex 1 for V /E R = {9.25, 11, 13}, where J is the 
tunneling matrix element. The rough constancy of this ra- 
tio accords with the general concept of a self-trapping thresh- 
old 03 1, but the absence of an initial internal expansion shows 
that the pure mean-field description is incomplete. 

The starting point for our mean-field theory is the Gross- 
Pitaevskii equation (GPE), 

/M f h 2 \ 

ih ~dt = \-m^ + v w+9m 2 J* ■ (i) 

Here = <b(x,t) is the matter wave field, g = i7 "^ , 
a s is the s-wave scattering length, V = Vo[sin 2 (7ra;/(i) + 
sin (pry/d)] is the optical lattice potential, and d is the lat- 
tice spacing. The interplay between free expansion along the 
tubes and the transverse lattice dynamics is accounted for self- 
consistently. 

The wave field (x, t) can be expanded in terms of Wannier 
functions along the optical lattice directions: 

V&t) = ^W n , m (x,y)<f> nm (z,t) . (2) 

n,m 

N n .m{t) = J dz\$> nm (z, t)\ 2 is the number of atoms 
at the site (n, m), satisfying the normalization condition 
Tim m Nnm = N, with N the total number of atoms. 
W n ,m (x, y) is the lowest band Wannier orbital localized at site 
(n, m). We assume, justified by the experiment, that only the 
first-band Wannier orbitals are occupied. We neglect any tem- 
poral dependence of the Wannier functions because the inter- 
well number/phase dynamics is much faster than the time as- 
sociated with the change in shape of such functions. 

Using the ansatz in Eq. [2] and integrating within each tube 
over the lattice directions (x and y), we get a set of coupled 
ID GP equations: 

( ti 2 d 2 2 \ 
\ 2M o 2 z J 

where $„ m = $ nm (z,t), J = f d 2 xW^ m [^j\7 2 - 
V]W n +i TO is the nearest neighbor tunneling, and U — 
I d- 2 x\Wn im \ 4 is the onsite interaction energy. We nu- 
merically solve the GPEs by discretizing them along z and 
calculating the kinetic energy term using a Fourier transform 
method. The widths obtained by assuming full adiabaticity 
during the lattice turn-on are larger than the experimentally 
observed widths. We therefore match the initial transverse 



width to the experimentally observed value, and adjust the ini- 
tial axial size so that the energy released in the axial expansion 
matches the experimental value. 




FIG. 3: (Color online) Transverse density distributions 12.4 ms after 
release for Vo = 13-Bjj. (a) r/=0. (b) rj=l. 

In the mean-field simulation [27] the central part of the 
atom cloud starts to expand at t = 0. Within several ms, 
self-trapping becomes evident with the formation of the pre- 
viously mentioned roughly square fort-like structure (see Fig. 
|3^). Outwardly moving atoms cannot pass the steep density 
gradient, and they reflect back. All the while, the self-trapped 
edges evolve, leading to a complicated and ever-changing den- 
sity distribution near the edges, and a slow spreading of the 
outer atoms. As the atoms expand axially along the tubes, 
the interaction energy decreases with time. Eventually the in- 
teraction energy differences between adjacent tubes become 
smaller than the tunneling energy, so that all self-trapping 
edges are lost and the cloud expands ballistically at a constant 
rate in the transverse direction. 

In order to compare the simulation with the experimentally 
measured RMS widths, we extract widths from the simulated 
distributions in a way that mimics the measurement process 
(i.e., we square the density distribution and convolve it with a 
Gaussian of RMS width 1.7 /im, corresponding to the imaging 
system resolution). This process removes fine features from 
the theoretical curves that are not resolvable in the experiment, 
but it does not guarantee that the theoretical distributions are 
well-fitted by Gaussians. In fact, the pure mean-field theory 
calculations retain clear evidence of self-trapped squares even 
after this processing [27]. The theoretical RMS widths are 
shown as blue dashed lines in Fig. [2] For all lattice depths pre- 
sented, the agreement is poor. For Vq = 7.25 Er, the mean- 
field calculation captures the qualitative dynamics, but over- 
estimates the expansion rate at early times. For Vq = 9.25, 11 
and 13 En, the calculation doesn't capture the qualitative dy- 
namics. It greatly overestimates the initial expansion rates, 
and predicts that the RMS width increases much more slowly 
during the later ballistic expansion than is observed. 

The mean-field treatment neglects thermal and quantum 
fluctuations, which can play large roles in the intermediate 
coupling limit. We note that it is hard to adiabatically trans- 
form a 3D gas into a ID gas in theory ll28l and, as described 
above, there is a small amount of irreversibility in turning on 
the lattice. Although the system is therefore not initially in its 
ground state, the observed self-trapping dynamics are insensi- 
tive to the details of the lattice turn-on and to the precise pu- 
rity of the initial BEC, making quantum fluctuations the bet- 
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ter candidate for the missing piece of the model. If we ignore 
tunneling and calculate the initial coupling strength, 7 11291 . 
of atoms near where atoms are predicted by mean-field the- 
ory to have a self-trapped edge, it varies from 0.2 to 0.31 for 
V = 7.25 E R and from 0.32 to 0.5 for 13 E R . Higher cou- 
pling implies larger phase variation along a tube. For small Vb, 
coherence among the tubes is maintained by tunneling. But as 
Vo is increased and the tubes become more ID, phase fluctua- 
tions along each tube will give rise to phase fluctuations from 
tube to tube. This scenario is reminiscent of the deconfine- 
ment transition that occurs when there is an additional axial 
ID lattice |3Q| . but the theory is less tractable out of equilib- 
rium and with no ID lattice, and as far as we know it has not 
been solved. 

In an attempt to incorporate quantum fluctuation we use 
an approximate TWA, which gives quantum fluctuations to 
leading order when quantum dynamics are expanded around 
the classical (Gross-Pitaevski) limit. To implement the TWA 
we need to solve the mean field equations supplemented by 
random initial conditions distributed according to the Wigner 
transform of the initial density matrix. Since the initial state 
is an interacting system, it is not simple to find the Wigner 
transform. Therefore, we accomplish an approximate TWA 
in an ad hoc way by introducing phase fluctuations into the 
initial conditions of our calculation by adding random phases 
among the tubes: (j> nm (t — 0) — > 2irrik, where n ,m(O) is the 
initial phase of the tube at lattice site (n, m), A is a uniformly 
distributed random variable between and 1, and 77 parame- 
terizes the strength of the phase fluctuations. Within each tube 
the gas is still described by mean-field theory. When 77 = 0, 
the initial conditions correspond to a fully coherent array of 
ID gases and when 77 = 1, we have a totally randomized, 
initially incoherent array. 

As long as 77 > 0.4, the squared and convolved theoreti- 
cal distributions fit fairly well to Gaussians. We average the 
results from 20 randomized phase implementations, and dis- 
play the RMS widths in Fig. [2] for various 77 values. For 
Vb = 7.25 Er, 77 ~ 0.4 best fits the early evolution, when 
there is self-trapping. For larger Vb, 77 = 1 clearly fits the early 
evolution best. Randomized phases cause initially random 
tunneling current directions throughout the tube array. Site- 
to-site density fluctuations rapidly develop, as seen in Fig.[3]5. 
Since there is no reflection from relatively sharp self-trapped 
edges, the acceleration of the cloud does not go negative (as 
it briefly does for 77=0), but instead asymptotically approaches 
zero from above. Self-trapping is still lost when the density 
drops below a critical value that depends on Vo. 

Although the modeled curves with the appropriate 77 are 
never too far from the observed ones for Vq = HE R and 
13-Er, the experiment shows a more sudden and more de- 
layed onset of ballistic expansion than the model, starting at 
a density that is at least three times smaller. We have per- 
formed calculations with fluctuations in the initial atom num- 
ber as well as phase in each tube and the curves are nearly 
unchanged. Assuming a normal error distribution, the p-value 
is only 2xl0~ 7 (6xl0~ 19 ) that the points between and 8 



(12) ms in the Vq — HEr (13Er) curves are consistent with 
mean field theory with maximal fluctuations. In contrast, the 
p-value is 0.45 (0.15) for them to lie on a horizontal line. The 
data suggests that the initial self-trapping is more complete 
than the model with maximal phase fluctuations predicts. 

Detailed two-body correlations within each tube may be re- 
sponsible. As the 2D lattice is turned on, correlations develop 
that make atoms avoid each other, at a kinetic energy cost, 
in order to avoid a larger mean-field energy cost. Since the 
wavefunction of a tunneling atom will not in general be ap- 
propriately correlated with the atoms in an adjacent tube, it 
would have to pay much of the mean-field energy cost that 
the correlations avoid. If that additional energy exceeds the 
tunneling energy, tunneling cannot proceed. This is a qual- 
itatively new self-trapping phenomenon which, in the strong 
coupling limit, is similar to Pauli blocking. It does not require 
density gradients, which is why the expansion of the central 
part of the tube bundle is fully suppressed until the overall 
density, and hence the mean-field cost of tunneling, gets suf- 
ficiently small. We speculate that such a process might tend 
to suppress mass transport whenever there is intermediate or 
strong coupling, regardless of the dimension lF3Tl[32l . 

In conclusion, we have measured the transverse expan- 
sion of quasi- ID arrays of atoms in the intermediate coupling 
regime. We find coarse agreement with a mean-field model 
when we introduce intertube phase fluctuations to the initial 
conditions. Modest phase fluctuations make the theory agree 
with experiment for Vb = 7.25 En for short times. For larger 
Vb, agreement with experiment is best with maximal initial 
phase fluctuations, which slow down the initial expansion. 
However, the experiment shows even less initial expansion, 
which suggests that there is strong self-trapping of a quali- 
tatively different type, perhaps due to microscopic quantum 
correlations. This behavior could in principle be captured 
by a full TWA calculation; however, more work is needed 
to quantitatively model tunneling behavior in this quantum 
many-body regime. 
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